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We investigated the finite temperature (T) phase transition for SU(A'^c) gauge theory with Nc = 4, 
6, 8 and 10 at lattice spacing, a, of 1/{6T) or less. We checked that these theories have first order 
transitions at such small a. In many cases we were able to find the critical couplings with precision 
as good as a few parts in 10''. We also investigated the use of two-loop renormalization group 
equations in extrapolating the lattice results to the continuum, thus fixing the temperature scale 
in units of the phase transition temperature, Tc- We found that when a < l/(8Tc) the two-loop 
extrapolation was accurate to about 1-2%. However, we found that trading Tc for the QCD scale, 
' ■'^Msi increases uncertainties significantly, to the level of about 5-10%. 
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I. INTRODUCTION 



Since the realization [Ij] that a non-trivial and tractable hmit is obtained for S\J{Nc) gauge theories when the gauge 
coupling, g, is taken to zero and the number of colours, is taken to infinity, keeping the combination g'^Nc fixed, 
there has been much work on this hmit . Most such work sums large classes of Feynman diagrams and therefore is 
^ . closely related to perturbation theory. The hope is that the hmiting theory and a small number of corrections in a 
' series in 1/Nc would allow us to understand the physically interesting theory with Nc ~ 3. Lattice calculations are 
, of help in testing this conjecture by making the connection from the other direction — by simulations and complete 
non-perturbative computations at finite Nc. They test whether a short series in 1/Nc for Nc> 3 extrapolates correctly 
I to the tractable limit of Nc — > oo. However, in order to test the continuum computations, one must also take the 
^ ' continuum limit of the lattice theories. This is the main thrust of this paper. 

T-H " The theory with Nc = 3 has been studied extensively before [3], and its continuum extrapolation using the renor- 
0^ malized weak coupling expansion has been studied and found to work Q . We shall have occasion to use these results 
at various points in this paper. The finite temperature transition has been studied before in 3-1-1 dimensions on 
, lattices with a = 1/(4T) for Nc — A These early studies found that the crossover from strong to weak coupling, 
0^ which is a lattice artifact, interfered with the finite temperature transition. Variant actions were invented to solve 
, this problem A modern solution which depends on today's vastly improved computational power is to just go to 
■ larger Nt with the simplest action. For larger Nc there have been some studies recently with Nt ^ 5, 6 and 8 For 
Nc — 4, A— has been extracted from data on the string tension in the Schrodinger functional scheme 

In this paper we investigate the continuum limit of the finite temperature deconfinement transition in SU(iVc) 
pure gauge theory for iVc > 3. The main thrust of our study is to control the approach to the continuum limit by 
performing simulations of the 3-1-1 dimensional theories at a succession of lattice spacings, a, and then using the weak 
coupling expansion for the extrapolation to zero lattice spacing. It turns out that with today's computational power 
it is quite possible to reach lattice spacings small enough for two-loop renormalization group equations (RGEs) to be 
useful for the continuum extrapolation. Indeed, at the lattice spacings that we use, even the one-loop flow is a good 
rough indicator of the continuum limit. 

In order to perform these precision tests of the continuum limit we performed lattice simulations of SU(4), SU(6), 
SU(8) and SU(IO) theories. In all cases we simulated theories with lattice cutoffs of a = 1/(6T) and l/(8r), and 
in some cases for even smaller lattice spacings, going down to lattice spacing of l/(12r) in one case. We performed 
finite size scaling studies, thus extrapolating to the thermodynamic limit of infinite spatial volumes, to check that 
the thermal phase transitions is actually of first order at lattice spacings a < 1/(6T). Coupled with the continuum 
extrapolations that we discuss next, this verifies earlier arguments about the order of the finite temperature phase 
transition in continuum theories with Nc > 3 ISl. 
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Through the finite size scaHng analysis we located the phase transition point with a statistical precision of a few 
parts in 10^. We found that the location of the phase transition point scales as expected in the limit of — > oo. 
With this precision we could test the two- loop RG flow to a statistical accuracy of a few parts in 10^. It turned out 
that at lattice spacing of a < l/(8Tc), the two-loop RGE is trustworthy in extrapolation towards the continuum, 
within 3ct of the statistical accuracy. In all this the quantity Tc is used to set the scale of measurements. 

Any test of a weak coupling expansion involves the choice of an RG scheme, i.e., a choice of a measurement used 
to define the running coupling in the gauge theory. If the perturbation theory is accurate, and all orders in the 
expansion are available, then the choice of the scheme is immaterial for any measurement. However, in all practical 
cases only a small number of terms in the weak coupling expansion are available. We found that for the determination 
of the temperature scale in terms of T^. the scheme dependence is statistically significant, but small in magnitude, 
being around 1-2%. This indicates that the lattice spacings used in our study are small enough for the use of the 
weak-coupling expansion. It seems likely that three-loop computations can improve matters. 

This could be the first indication that non-perturbative lattice computations for Nc > 3 are at a point where they 
are more reliable than the perturbative series needed for the continuum extrapolation. Needless to say, one could 
just push the non-perturbative lattice simulations to smaller and smaller a until the running coupling (at the scale of 
a) decreases significantly and the available perturbative expansions begins to be more accurate. However, it is more 
cost-effective to develop the perturbation theory to higher order. 

In performing a weak coupling expansion the scale of choice is one which defines how fast the coupling changes 
asymptotically when measured at two different length scales. This intrinsic scale of QCD is called Ajj^. We found 
that the determination of Aj^ in terms of the non-perturbatively determined scale Tc is quite uncertain. While the 
statistical errors are under control, the scheme dependence is quite large. Our observations seem to indicate that one 
needs smaller lattice spacings to stabilize the transformation from Tf. to A-jjg. 

This paper is structured as follows: in the next section we discuss the technicalities of the lattice simulations. 
Following this we present our results for the finite temperature transition and its extrapolation to the thermodynamic 
limit. Next, we discuss the continuum limit, the setting of the temperature scale and the extraction of Aj^. The 
final section contains a summary of our results. Some parts of our results have been reported earlier in conference 
proceedings ^10.] . 

II. SIMULATIONS, MEASUREMENTS AND OTHER TECHNICALITIES 



Nt 


iVc = 4 Nc = G Nc = S Nc^ 10 


4 


12, 16, 18, 20, 24 12, 16, 20 


6 


16, 18, 20, 22, 24 14, 16, 18, 20, 24 16* 16* 


8 


22, 24, 28, 30 20, 24* 16* 16* 


10 


24 24* 


12 


24 



TABLE L For each Nc and Nt the values of Ns used in the simulations are given. Runs which are exploratory are marked by 
asterisks. The remaining runs are meant to yield precision data; for these the details of the statistics are given in Table IVll 
Zero temperature runs were performed for Ns — Nt — 16 for all Nc and Ns — Nt = 24 for A'^c = 4 and 6. 

In this study we use the Wilson action — 

S = (3J2 [l-ReP^,(z)], (1) 

where -P^;y(i) is the trace of the product of SU(A'c) valued link matrices, U, around a plaquette, starting from the 
site i and touching the site i + + v. The trace is normalized by a factor of Nc, so that by this definition the trace 
of an unit matrix is unity. The lattices have size Nt x in units of the lattice spacing, a. The physical extent of 
the lattice is aNt = 1/T and £ = aNg = C,/T where C — Ng/Nt is called the aspect ratio. Increasing ( at fixed T 
corresponds to increasing the volume, V — The bare gauge coupling is = 2Nc/P. 
The partition function. 



Z{V,T) = / VUe-^^^\ 
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FIG. 1: The average plaquette as a function of the bare coupUng for TVc = 4 on a 16* lattice. Above /3 = 10 the strong coupling 
series no longer predicts (P) accurately, and the theory crosses over to the weakly coupled phase. 
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FIG. 2: The average plaquette as a function of the bare coupling for Nc = 4 and 6 on various lattice sizes. For SU(4) on a 
lattice with Nt — 4: there is a jump in (P) at /3c ~ 10.48 where the first order thermal phase transition occurs with a jump in 
L. However, at larger Nt there is no jump in (P). For SU(6) there is a jump in (P) at all Nt. For A^t = 4 the jump occurs at 
the thermal phase transition, but at all other Nt the bulk and thermal transitions are decoupled. 



is sampled using a Monte Carlo procedure in which over-relaxation steps are mixed with heat-bath updates. A large 
fraction of the CPU time is taken up in the computation of the product of matrices connecting to a given matrix 
(called staples). This computation scales as N^, since the time is dominated by the multiplication of Nc x Nc matrices. 
Therefore, for each computation of a staple, it would make sense to update each of the Nc{Nc — l)/2 SU(2) subgroups 
of SU(A^c) a fixed number of times [llj. When we update all SU(2) subgroups once in every step of a composite sweep 
which contains three steps of over-relaxation per step of heat-bath, then about 50% of the CPU time is spent in the 
computation of staples, about 33% in the over-relaxation update, and about 12% in the heat-bath. The rest of the 
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time is spent in the measurement of plaquettes and Polyakov loops. These fractions are almost independent of Nc, 
whereas the actual CPU time per link update scales very close to N^. It was argued earlier [l^ that in an optimum 
hybrid over-relaxation algorithm the number of over-relaxation steps should be increased linearly with Ng. If this 
were to be done, then relatively less time would be spent in computing staples, resulting in more optimal use of CPU 
time. 

We performed simulations of four theories. An overview of the runs is given in Table [J and its caption. Almost 
all zero temperature runs collected statistics of several tens of thousands of composite sweeps, and most runs have 
statistics of over half a million composite sweeps. The statistics of a set of measurements should actually be judged by 
the auto-correlation time, Ti„t, since the error in a measurement, E, is related to the variance of the measurements, 
(7^, through the formula E'^ = Tintcr^ /N where N is the number of measurements. Auto-correlation functions of the 
plaquette at T = show that Tint varies between approximately 1 and 10 sweeps. Since we study first order phase 
transitions. Tint in the transition region for the order parameter, L, is closely related to the number of tunnelings 
between different phases [l^ . The statistics collected close to the transition region are summarized in Table IVII 

SU(A^c) theories with the action in eq. ([T]) exhibit a bulk transition when Nc is large enough. This transition can 
be monitored in zero temperature simulations using the expectation value of the plaquette, (P), where 

and d = 4 for our purposes. On the small-/? side of the transition, one expects the strong coupling series for (P) 
to work; this is an expansion of (P) in powers of 0^ [l4l|. At larger (3 one expects renormalization group running 
of (P) [15[. For Nf, — 4 the change from strong to weak coupling behaviour is fairly smooth, with a cross over in 
the vicinity of /3 = 10.2 (see Figure [T]). The strong coupling side has little to do with continuum physics. We study 
thermal physics on the weak-coupling side of this crossover, where, as we show in Section HVl the continuum limit can 
be taken. 

The largest finite volume effect at T = is expected to occur when the lattice sizes are such that a spurious 
deconfinement transition takes place [IB]. At any given bare coupling there is a critical such that for N'^ 

lattices with N > N^{/3) one expects small finite size effects. These small effects are expected to scale as exp(— ^toq) 
where toq is the mass of the lowest glueball. For SU(3) pure gauge theory, this mass is very high compared to the 
deconfinement temperature Tc- If this happens also for Nc > 3, then one expects that finite size effects should be 
smaller than of order ex-p[~N/N^{i3)]. That finite volume effects are indeed small at T = is borne out by the data 
in Table Eni 

The finite temperature transition was monitored using the order parameter Polyakov loop, L, where 

^ = ]^ETrn%, (4) 

I t=l 




10.4 10.6 10 



FIG. 3: L and (P) at functions of /3 on 6 x 18^ and 6 x 24^ lattices for SU(4). A rapid change at = 10.78 is seen in L, 
whereas (P) is continuous. For A*'t > 4 the bulk and thermal transitions are decoupled for all Nc, as shown by this kind of 
observation. 
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FIG. 4: A first order transition, i.e., tlie coexistence of phases, witli different vaiues of L, is signaled by a multi-peaked histogram 
of |L| and the fact that the scatter plot of L in the complex plane shows 5 well developed dense regions — L = Q and four 
complex values of L. Here we show these features at a coupling where all five coexisting phases for the SU(4) theory have large 
weight. 
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FIG. 5: Re-weighting analyses of X-c for SU(4) gauge theory shows that on the larger lattices the maximum scales with the 
lattice volume, , indicating a first order phase transition. The analysis is shown for both Nt = 6 and 8. 

where the sum over sites, i, is restricted to all spatial sites. The order parameter jumps from a zero value at small 
temperature to a finite value at the thermal transition, signaling deconfinement. The thermal transition is of first 
order in all the simulations presented here. We found that for SU(4) and SU(6) gauge theories the finite temperature 
transition and the bulk transition interfere for iVt = 4 (see Figure [2]) . This is a known phenomenon 0, 0] . Since the 
bulk transition must occur at a fixed lattice spacing, it is natural to expect that by changing Nt the bulk and the 
thermal transitions can be decoupled. It was found Q that at larger Nt these transitions do separate out (see Figure 
|3]). Therefore, our strategy in this paper is to study larger iVt, where the thermal transition is in the weak coupling 
regime, and use these studies to take the continuum limit. 



III. THE DECONFINEMENT TRANSITION 



The abrupt change of L shown in Figure [3] indicates that the finite temperature transition could be of first order. 
Clear evidence of the coexistence of phases labeled by the value of (L) is obtained from the distribution of i. In 
simulations of the SU(4) theory close to (3c we found that the system is equally likely to be in the phase with L = and 
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FIG. 6: Finite size scaling for SU(4) gauge theory for Nt = 6 (boxes) and 8 (circles). The first panel shows the maximum of 
Xi/(Vr^) as a function of 1/(VT''). The second panel shows /3c as a function of 1/{VT'^). On the largest spatial volumes, 
the maximum scales as V , as expected for a first order phase transition. On the same volumes l3c reaches a limit which is its 
thermodynamic value. 




FIG. 7: Re-weighting analyses of Xl for SU(6) gauge theory shows that on the larger lattices the maximum scales with the 
lattice volume, V, indicating a first order phase transition. Also shown is a multi-histogram analysis on the largest lattice with 
seven and three input simulations, demonstrating the stability of the estimate of /3c- 



in four phases with the same (|L|) but different phase angles (Figure^]). Hence the histogram of \L\ shows two peaks, 
one close to zero and another elsewhere. A scatter plot of L measured on each gauge field configuration also shows 
four distinct populations. All these observations are consistent with a first order phase transition. The extraction of 
the jump in (L) at Tc needs the renormalized Polyakov loop ^Y?\ and hence lies beyond the scope of this study. 

For more accurate determination of Pc we defined this coupling by the position of the maximum of the susceptibility 
of \L\- 

X. = N! {{\L\') - {\L\')} . (5) 

For the exploratory runs marked in Table |T1 Pc is estimated from the position of the peak of the values of Xl found 
in a scan over /3, and its quoted error is the spacing in the scan of (3. In all the remai ning cases, the objective was 
precision, and maximum of Xl was determined through multi- histogram re- weighting jlSf . The errors on Xl were 
defined through a bootstrap procedure combined with the re-weighting. Such analysis requires very large statistics, 
which is available to us, as shown in Table IVTl 

A final verification of the order of the transition and the determination of the critical coupling, (3c, require finite 
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iVt = 4 iVt = 6 Nt^8 Nt = 10 Nt = 12 



3 



5.6925 (2) 5.8940 (5) 6.0609 (9) 



10 



4 



6 



8 



10.788 (1) 11.078 (1) 11.339 (4) 11.552 (17) 
24.838 (1) 25.470 (3) 26.0 (1*) 
44.7 (2*) 45.8 (2*) 

70.5 (15*) 73 (2*) 



TABLE II: The critical couplings, /3c, for the first order thermal phase transition for different Nc and different temporal lattice 
sizes, Nt. Error estimates which are marked by an asterisk are not statistical, as discussed in the text. The results for Nc — 3 
were found in [3]. For Nc > 3 the transition for Nt — 4 falls in the region of the strong to weak coupling cross over, making it 
hard to distinguish the bulk from the thermal phase transition. 



when is large enough. Also, the position of the peak, which is our estimate of /3c at finite volume, should scale as 



as one approaches the thermodynamic limit, Ng — s- oo. A different definition of /3c (-^s), such as the one where the 
A'c + 1 different peaks in L have equal weight, could give a different result at finite through a change in 6. Finite 
volume scalings as in eqs. (0 [7]) were observed in SU(3) gauge theory [^. The asymptotic region sets in when the 
lattice size is much larger than the longest correlation length in the system. In this asymptotic region one expects 
exponentially slow sampling through a standard Monte Carlo procedure. Tint oc exp{aV^^^) |2^. As a result, one 
might expect that as the transition becomes stronger it becomes harder to do a finite size scaling analysis because of 
an increase in a, but the asymptotic finite volume corrections, 7 and S, also become relatively smaller. 

The variation of Xl with /? obtained through a bootstrap multi-histogram analysis is shown for the SU(4) theory 
with Nt = 6 and 8 in Figure In each case, it turns out that on the two largest lattice used, i.e., those with the 
aspect ratio, C > 3, one enters the region of asymptotic finite size scaling where the formulae in eqs. ([H [7]) become 
applicable. A more detailed view of the peaks of the susceptibilities, shown in Figure [51 bears this out. 

From this analysis, we can estimate the position of the peak of Xlj i-S-, /3c) since this is very stable on the largest 
lattices used. This stability is apparent in Figure [SI Previous measurements of S for A^t = 5 when extrapolated to 
larger Nt using the formulae used there predict much larger shifts than we observe. Our results for f]c are collected 
in Table mi 

SU(6) follows the same trend. For all Nt, one has all the qualitative features of a strong first order phase transition — 
multiple coexisting phases (6 ordered phases and one disordered in this case) and long auto-correlation times deter- 
mined by the tunneling rate from one phase to another, growing rapidly with volume. The phase transition is even 
stronger than SU(4), and a finite size scaling analysis is more delicate. 

In Figure [3 we show the multi-histogram reweighting analysis for SU(6). Note that the aspect ratios used in this 
analysis are smaller than those for SU(4). This is forced on us because the transition is stronger, and therefore Tint 
grows faster with V. In fact, statistical problems already begin to show up at the largest V for Nt = 6; the run at 
P = 24.85 has statistically too few tunnelings, since it lies right at the edge of the region of metastability for these 
lattices. For this system we examined the stability of the analysis through the comparison of the multi-histogram 
method with seven and three histograms. As shown in Figure [3 the peak is unambiguously determined, since the 
values of xT the two analyses are compatible, as are the estimates of /3c. The reason for the absence of a large 
systematic error at the peak is that the scan in /3 is fine enough, so that there are enough other histograms to 
compensate for the one which is improperly sampled. 

Our simulations of SU(8) and SU(IO) pure gauge theory at finite temperature were purely exploratory, being 
restricted to a single volume at each Nt. The value of (3c that we estimate, along with the error bounds given by the 
scan in (3 are quoted in Table [Til 

In Figure [51 we plot these results as a function for A'^c at fixed lattice spacing a — l/{NtTc) for A't = 6 and 8. We 
see that a good description of our observations is obtained by a two-term extrapolation to the large- A'c limit — 




(6) 



Pc{nI)^Pc^5n-^ + o{n;^), 



(7) 



N^c ^* Nr 



(8) 



8 




FIG. 8: /3c for different number of colors, on A't = 6 and 8 lattices. 



The quantity /3, is expected to increase without bound as Nc — > oo. The first correction term, of order provides 
a sufficient description of the data even at Nc = 3. This scaling check shows that for each cutoff, a — l/{NtT) one 
has a large Nc theory which is non-trivial in the limit g^Nc fixed, i.e., /3/Nc fixed. 

Note that in the best cases we have achieved accuracies of a few parts in 10'' in the measurement of /3c. Next we 
turn to the continuum extrapolation of these measurements and the determination of the temperature scale. 



IV. RENORMALIZED COUPLING AND THE TEMPERATURE SCALE 



Pure gauge SU(-/Vc) theory contains a single dimensionless parameter, the coupling, as — /At:. Quantum cor- 
rections change this into a scale. This can be specified explicitly, as the parameter A, or implicitly, as the value of 
the running (renormalized) coupling as{fJ-) at a chosen momentum scale /i. At scales where as is small, perturbation 
theory is expected to work. In that case, changes of the scale of measurements can be accomplished through the use 
of perturbation theory. In particular, extrapolation of results to the continuum can then be done with ease. 

The two- loop RGE can be integrated to trade the running coupling as(M) for a mass scale, 

aK = kR\ — I I, where exp(-.x/2)a;^i/(2''"), (9) 

where k depends on the coupling as that enters into these equations. This coupling is measured by some operator 
dominated by the ultraviolet scale 1/a. Each such definition of as defines an RG scheme. The function R is obtained 
by integrating the two-loop beta function, 

^(3)=M^ = -/5o.9'-/?i5', (10) 

where /3o and /3i are well-known [l^]. These coefficients are independent of the scheme. Since T = l/{aNt), and we 
have a determination of Tc for different iVj, by making appropriate lattice measurements of as we can measure the 
temperature scale, T /Tc. At the same time, one could use eq. © to determine the QCD scale A-jjs" in terms of Tc. 

In order to complete this process, we need to define as- Two schemes are easily implemented on the lattice. One is 
the V scheme , in which the potential extracted from Polyakov loop correlations is used to define the renormalized 
coupling. Equivalently, at two-loop order accuracy, the weak-coupling expansion of the plaquette [23j can be inverted 
to find Uv — 



ln(P) =TTCFa^{q) 



\2tt \ aq 



In ay{q) 



(11) 



where Cp = (N^ — l)/(2-/Vc) and q = k/a, where k is the same number which is used in eq. ©. In this scheme 
k = 3.4018 [15|. Since (P) is easily measured and needed for thermodynamic quantities, we prefer to use eq. (jll[) as a 
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definition of av rather than through a separate measurement of the potential. The other definition is the E-scheme, 
in wliicli the couphng is defined from the plaquette through the formula 



1- (P> =7TCFaE{q), 



(12) 



where q = 1/a, i.e., k = 1. If the weak coupling expansion were exact, and known to all orders, then there would be no 
difference between the couplings determined in these two schemes at any cutoff, provided that ay (or a^) were small 
enough. Since this is not the case, one must explore RG scheme dependence. A third scheme that we utilize is the 
MS scheme defined through dimensional regularization of the continuum perturbation theory. The known expansion 
of Uv in terms of ajf^ [24] is used to obtain the latter using the two-loop relation 



1 



2N, 



(13) 



where q' — exp(— 5/6)g In other words, k — 1.4784 for the MS scheme. 

The values of the plaquette at zero temperature are measured on the grid of (3 shown in Tables IVIII and IVIIII They 
are obtained at other points using Lagrange interpolation with polynomials of orders between 1 and 4, and through 
a cubic spline interpolation. By using such a variety of interpolation schemes we quantify the systematic error in 
the interpolation at any fi as the widest dispersion between these schemes. For SU(4) and SU(6) on lattices with 
Nt > 6, this systematic error is smaller than, or of the same order as, the statistical error in the measurement of the 
plaquette. For SU(8) and SU(IO), the systematic error is larger than the statistical error. These lead to statistical 
and systematic errors in the determination of the running coupling of the order of a few parts in 10^. However, when 
we determine a scale, the largest error is that which comes from the determination of /3c- 

A test of the weak coupling expansion for the scale, and the scheme dependence in this is provided by using the 
determination of /3c for one Nt to predict that at a different Nt- Since we have measurements for A'^f = 6, 8 and 10 
for SU(4) and SU(6), we have chosen to examine the temperature predicted by the one-loop and two- loop RGEs for 
the A't = 6 and 10 lattices at the /3c corresponding to the Nt — 8 lattice. This is shown in Table IIIII Note that the 
error of roughly one part in 10** in the determination of /3c translates into an error of about one part in 10'^ in the 
determination of the temperature scale in the range of temperatures we explore here. Since the accuracy of this error 
estimate is important in our later reasoning, we performed it by two different methods: first by the usual methods of 
propagating errors, and then again through a bootstrap analysis. The two errors agreed within 10%, indicating that 
the estimates are robust. The errors quoted in Table [TTTl come from the bootstrap analysis. 





Nt 


Scheme 


2-loop 


1-loop 


4 


6 


E 
V 
MS 


1.29709 (167) (7) (1) 
1.30782 (174) (8) (1) 
1.30057 (169) (7) (1) 


1.32333 (184) (8) (1) 
1.35339 (204) (9) (2) 
1.35068 (202) (9) (1) 




10 


E 
V 
MS 


0.80885 (265) (5) (0) 
0.80442 (270) (5) (0) 
0.80757 (267) (5) (1) 


0.79632 (280) (5) (0) 
0.78381 (294) (5) (1) 
0.78517 (292) (5) (0) 


6 


6 


E 
V 
MS 


1.30504 (165) (10) (3) 
1.31675 (172) (10) (4) 
1.30916 (167) (10) (3) 


1.33190 (181) (11) (4) 
1.36476 (200) (12) (4) 
1.36130 (198) (12) (4) 




10 


E 
V 
MS 


0.81884 (2993) (4) (1) 
0.81435 (3052) (4) (1) 
0.81740 (3011) (4) (1) 


0.80695 (3161) (4) (1) 
0.79452 (3324) (4) (1) 
0.79583 (3307) (4) (1) 



TABLE III: The values of T/Tc at the coupling /3c{Nt = 8), for SU(7Vc) gauge theory for different At, in different RG schemes, 
and at different loop orders. The entries give the central value, the statistical error propagated from the uncertainty in l3c{Nt), 
the statistical error from plaquette measurements, and systematic errors from interpolations of plaquette values. For At = 6 
the exact non-perturbative result is T/Tc = 1.33, and for At = 10 it is T/Tc = 0.80. 

Some systematics visible in Table IIIII is worth explicit comment. The one-loop RG already is close to the exact 
result, but in all cases performs worse than the two-loop RG. This is expected. Also the RG flow between f3c{Nt = 10) 
and (3c{Nt = 8) is better than that between f3c{Nt = 6) and PdN — 8). Indeed, for SU(4), where the test is most 
stringent, the former agrees with the exact non-perturbative result to about 1.5a in the V-scheme, and to about 3-4fT 
in the other schemes. The implication is that A^t = 8 is already in a regime where the weak coupling extrapolation to 
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CO 

II 



^ 1.4 




FIG. 9: One-loop renormalization group flow for Nc — 4: and 6 in the V-scheme. The data on /3c for A^t = 5 is taken from [t!]. 
If the RG were adequate, then the curves for different Nt would lie on top of each other. The accuracy of the one-loop flow 
improves with increasing Nt- 




FIG. 10: Two-loop renormalization group flow for Nc = 4 and 6 in the V-scheme. The data on /3c for A^'t = 5 is taken from 0]. 
Good scaling behaviour is obtained when the curves for different Nt lie on top of each other. Excellent scaling is obtained for 
iVi = 8 and 10. 



the continuum works, but Nt — 6 may be just a little outside this regime. The scheme dependence is most significant 
when the RGE is least reliable, but is less than 1% in all cases. 

The one-loop temperature scale is shown in Figure [5] in the V-scheme for a large range of lattice spacings. While 
this works reasonably well, the improvement in going to two-loops, shown in Figure [TOl is obvious. The cutoff effects 
are small on this scale already for Nt = 5. The two-loop temperature scale in the V-scheme for Nt = 6 and 8, are 
collected together in Table IIXI for future reference. As discussed already, the scale for A't = 8 is more reliable, and 
should be used for extrapolations. The scale for Nt — 6 serves to give a rough indication of the kind of systematic 
errors to be expected: as one can see the difference between these two scales is roughly of 2%. 

We have argued here that the continuum extrapolation from Nt = 8 or 10 can be performed using the weak coupling 
expansion. In weak coupling the reference scale that is used is A— and not Tc. Our analysis above gave us the scales 
Ae and Ay. These can be converted into A^ [26l. [27|. 



cxp 



1.622268-^) 
11 V 4A^cV 



^7 



exp 



31 
66 



(14) 



Using this and the determinations of /3c, we can convert the non-perturbative scale Tc into a specification of A-jj^- in 
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iVc 


E- scheme 


V-scheme 
Tc/A— C2 


MS-scheme 
Tc/A— C2 


3 
4 
6 


1.16 (2) 2.7 (7) 
1.198 (1) 2.49 (5) 
1.193 (1) 1.98 (9) 


1.11 (2) 1.9 (7) 
1.129 (2) 1.61 (5) 
1.120 (2) 1.05 (9) 


1.17 (2) 2.5 (7) 
1.203 (2) 2.23 (5) 
1.193 (2) 1.67 (9) 



TABLE IV: Fitted parameters for the non-perturbative beta function in the form of eq. (|15|l . Here T^/ A-^s must be considered 
as a formal fit parameter. Data from all available lattice spacings l/(12Tc) < a < l/(6rc) have been used. For TVc = 3 data 
from a = l/(4rc) has also been used. 





E-scheme V-scheme MS-scheme 


3 


1.19 (3) 


1.12 (3) 


1.20 (2) 


4 


1.235(1) 


1.153(1) 


1.236(1) 


6 


1.222(1) 


1.135(1) 


1.217(1) 


8 


1.26(6) 


1.17(5) 


1.25(6) 


10 


1.48(41) 


1.38(39) 


1.48(41) 


oo 


1.22 


1.13 


1.22 



TABLE V: TcJ Ajfg in the contirmum limit of SU(A'^c) gauge theory for Nc ~ 3,4,6,8, in different schemes using two-loop RGE 
for a < l/{STc). These values of Tc/A-j^ are appropriate for use in a two-loop computation. 

three different schemes. 

The test of two-loop RGE in Table IIIII showed that this was fairly accurate already at the lattice spacing corre- 
sponding to Pc for Nt = 8 and 10. Therefore it is no surprise that to the same degree of accuracy the ratio Tc/A— is 
constant when evaluated at Nt > 8. However, the scheme dependence is much larger for A-jj^ than for the temperature 
scale. This happens because the renormalized coupling is not small enough for the eqns. (fT4|) to hold. One could 
correct these formulae by explicitly including two-loop or higher order correction terms. However, then the ratio 
Tc/Ajfs would depend on the scale a. This can be avoided only when as becomes substantially smaller. However, 
since as runs logarithmically with a, that would imply that one has to use lattice spacings which arc about 10 times 
smaller. This is currently outside the reach of our computational abilities. 

For SU(3) the range of accuracy of the RGE can be extended by including into it corrections of order ^22,]. This 
seems to be possible for iVc > 3 too. When we compare Tc/Ajf^ extracted for all Nt, it seems possible to fit this to a 
simple variation. In principle this term can be used to add an ©(a^) correction to the two-loop beta function by 
changing R{x) in eq. ([9]) to R{x)[l + rj/Nt]. We evaluate these corrections by a fit to the lattice spacing dependence 
of Tc/A— which renders this ratio flat in the whole range l/(10Tc) < a < l/{5Tc}, i.e., we choose the fit form 

Statistically significant results can only be obtained for Nc < 6. Our results for the fit are given in Table ITVl One 
can compare these with the estimate Tc/ Ajts = 1.187 ± 0.009 obtained by combining estimates of Tc/ \fa (where cr is 
the string tension) and y/a / k-^ reported in 0, [l^. The estimation of \fo / k-j^ removes 0(c?) corrections, as we do 
here. We note that such a term sums many different types of corrections and amounts to a phenomenological fit of 
the beta function, i.e., gives what is called the non-perturbative beta-function. For this reason it cannot be regarded 
as a test of scaling. 

The continuum values for Tc/kjrsi obtained assuming that this ratio is constant for lattice cutoffs a < l/(8Tc) are 
collected in Table [V] Note that there are large and statistically significant differences between these results and those 
in Table IIVI Since the latter results constitute a check of two-loop RGE, and the best possible extraction of A-j^g, 
they are to be preferred for this purpose. For SU(3) we have performed a re-analysis of the data which was used in 
Q without the 0{a?) terms from [2^. This makes the analysis uniform for all Nc- Note that the dependence on Nc 
is weak. We have added indicative values of this ratio extrapolated to the limit Nc ^ oo. Since a statistical analysis 
is not possible, we have not added error bars to this extrapolation. Note that the strong scheme dependence, which 
we discussed before, propagates to the Nc ^ oo hmit. 

In summary, the two-loop renormalization group equations work well for a < l/(8Tc), i.e., at the level of 1-2%. 
Since the largest part of this uncertainty stems from the RG scheme dependence, higher order corrections in the 
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perturbation series for the plaquette could easily improve this description. However, trading the non-perturbative 
scale Tc for the perturbatively determined scale Kjj^ is not yet possible to better than 5-10%. Improving this would 
require using lattice spacings which are beyond reach today. 

V. CONCLUSIONS 

In this paper we studied the finite temperature phase transition in SU(4), SU(6), SU(8) and SU(IO) pure gauge 
theories at several lattice spacings and extrapolated the results to the continuum. In all these theories at large 
lattice spacing, a ~ l/(4Tc), a lattice artifact called the bulk phase transition prevents a simple study of finite 
temperature physics. The order parameter of the bulk transition is the plaquette average, (P), whereas that of the 
finite temperature transition is the Polyakov loop expectation value, L. The bulk transition is expected to occur 
at a (approximately) fixed lattice spacing. We studied these theories at smaller lattice spacings, a < l/(6Tc), and 
found that in all cases the finite temperature phase transition can be studied without any interference from the bulk 
transition (see Figure [U for example). More details are reported in Section [TTl 

We found a first order finite temperature transition for all these theories. This was established not only by clear 
signals of multiple coexisting phases labeled by different values of L, but, in several cases, also by finite size scaling 
tests. These studies and also multi-histogram reweighting at fixed volumes allowed us to locate the phase transition 
with precision which was in many cases as good as a few parts in 10^. Our results on the finite temperature phase 
transition are given in Section IIIIl and the locations of the phase transition are collected together in Table HIl 

We investigated the continuum extrapolation of our lattice results and found that when the lattice spacing is 
a < l/(87c) then the two-loop RGE can be used to take the continuum limit. In order to do this one has to use 
a definition of the renormalized (running) coupling, called an RG scheme. We found that when the location of the 
phase transition at one lattice spacing is used to predict that at another, then the dependence on the RG scheme is 
small (see Table llll| : the statistical precision is about one part in 10"^, but the scheme dependence is about 2%. This 
allows us to construct a temperature scale with this degree of precision using the non-perturbatively obtained mass 
scale, Tc- Since the scheme dependence is the largest part of the uncertainty, higher order corrections will reduce this 
error. This is the first instance of a large Nc lattice calculation which has reached precisions good enough to test the 
state of the art in the weak coupling expansion. Details of these tests can be found in Section IIVI One useful result 
is the determination of the temperature scale in SU(4) and SU(6) gauge theories (Table llX|) . 

We tried to use two-loop perturbation theory to trade the scale Tc for the scale A-jj^ which is more commonly used 
in weak coupling expansions, and found that the scheme dependence becomes significantly more pronounced. The 
extraction of A-jt^- by this means gave statistical errors comparable to the temperature scale, but RG scheme depen- 
dence of about 10%. A large scheme dependence in trading a non-perturbative scale such as Tc for the perturbative 
scale Ajj^ is bound to persist in all foreseeable lattice computations. 

We found that two results can be easily extrapolated to the limit Nc ^ oo. The location of the critical point at 
fixed lattice spacing a — l/{NtTc) goes as /3, -I- ©(l/iV^) for Nc > 3 (see Figure [8|). For Tc/ A-jjs the series could 
be shorter; we find no statistically significant dependence of Tc/ A— on Nc in any of the three RG schemes that we 
studied (see Table |V|. 

These computations were carried out on the CRAY-Xl of the ILGTI in TIFR, and on the workstation farm of the 
Department of Theoretical Physics, TIFR. We would like to thank Ajay Salve for technical support. 



13 



[1] G. 'tHooft, JVucJ. Phys., B 72 (1974) 461. 

[2] See, for example, E. Brezin, C. Itzykson, G. Parisi and J.-B. Zuber, Comm. Math. Phys., 59 (1978) 35; 
E. Witten, JVucJ. Phys., B 156 (1979) 269; 
T. Eguchi and H. Kawai, Phys. Rev. Lett., 48 (1982) 1063; 

S. Coleman, in Aspects of Symmetry, Cambridge University Press, 1985, Cambridge, UK; 

E. Brezin and S. Wadia, The large N expansion in quantum field theory and statistical physics. World Scientific, 1993, 
Singapore; 

M. J. Teper, Phys. Rev., D 59 (1999) 014512. 
[3] Y. Iwasaki et al, Phys. Rev. Lett., 67 (1991) 3343; 

G. Boyd et al, NucL Phys., B 469 (1996) 419. 
[4] S. Gupta, Phys. Rev., D 64 (2001) 034507. 
[5] A. Gocksch and M. Okawa, Phys. Rev. Lett., 52 (1984) 1751; 

G. G. Batrouni and B. Svetitsky, Phys. Rev. Lett, 52 (1984) 2205; 
M. Wingate and S. Ohta, Phys. Rev., D 63 (2001) 094502; 
R. V. Gavai, JVucJ. Phys., B 633 (2002) 127. 
[6] S. Datta and R. V. Gavai, Phys. Rev., D 62 (2000) 054512. 
[7] B. Lucini, M. Teper and U. Wenger, Phys. Lett., B 545 (2002) 197; 
B. Lucini, M. Teper and U. Wenger, J. H. E. P., 0401 (2004) 061; 
B. Lucini and M. Teper, J. H. E. P., 0502 (2005) 033. 
[8] B. Lucini and G. Moraitis, Phys. Lett., B 668 (2008) 226. 
[9] B. Svetitsky and L. G. Yaffe, JVucJ. Phys., B 210 (1982) 423; 
B. Svetitsky, Phys. Rep., 132 (1986) 1. 

[10] S. Datta and S. Gupta, arXiv: 0906. 3929] 

[11] Ph. de Forcrand and O. Jahn, e-print h ep-lat/0503041 

[12] U. Wolff, Phys. Lett., B 288 (1992) 166. 

[13] A. Billoire et al., JVucJ. Phys., B 358 (1991) 231. 

[14] J. M. Drouffe and K. J. M. Moriarty, Phys. Lett., 108 B (1982) 333. 

[15] G. P. Lepage and P. B. Mackenzie, Phys. Rev., D 48 (1993) 2250. 

[16] S. Datta and S. Gupta, Phys. Lett., B 471 (2000) 382. 

[17] S. Gupta, K. Huebner and O. Kaczmarek, Phys. Rev., D 77 (2008) 034503. 
[18] A. M. Ferrenberg and R. H. Swendsen, Phys. Rev. Lett., 63 (1989) 1195. 
[19] C. Borgs and R. Kotecky, Phys. Rev. Lett., 68 (1992) 1734; 

S. Gupta, A. Irback, M. Ohlsson, JVucJ. Phys., B 409 (1993) 663; 

A. Billoire, NucJ. Phys. Proc. SuppL, 42 (1995) 21. 
[20] S. Gupta, Phys. Lett., B 325 (1994) 418. 
[21] B. A. Berg, Int. J. Mod. Phys., C 3 (1992) 1083. 
[22] C. Amsler et al, Phys. Lett., B 667 (2008) 1. 
[23] T. Klassen, Phys. Rev., D 51 (1995) 5130. 
[24] M. Peter, JVucJ. Phys., B 501 (1997) 471. 

[25] S. J. Brodksy, G. P. Lepage and P. B. Mackenzie, Phys. Rev., D 28 (1983) 228. 
[26] B. Alles, A. Feo and H. Panagopoulos, JVucJ. Phys., B 491 (1997) 498. 
[27] M. Luscher and P. Weisz, JVucJ. Phys., B 452 (1995) 234. 
[28] C. AUton, liep-lat/9610016 

R. G. Edwards, U. Heller, T. Klassen, NucJ. Phys., B 517 (1998) 377. 
[29] C. AUton, M. Teper and A. Trivini, J. H. E. P., 0807 (2008) 021. 



14 



APPENDIX A: SOME DETAILS 

Some details of the simulations and detailed tables of some of our results are collected in this appendix. 











= 4 












= 6 




Nt=6 


iVt =8 




iVt =8 
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Statistics 


Tint 
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Statistics 


Tint 
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Tint 
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Statistics 


Tint 


14 














24.80 


2.14 


7354 






















24.81 


1.94 


9391 






















24.82 


2.05 


10538 






















24.83 


1.94 


12349 






















24.84 


1.94 


12113 






















24.85 


2.05 


9428 






















24.86 


2.01 


7747 








16 


10.77 


1.56 


2367 








24.80 


2.26 


4388 










10.78 


4.5 


3188 








24.81 


3.88 


12602 










10.79 


1.56 


3780 








24.82 


2.29 


14827 










10.80 


2.6 


2566 








24.83 


2.74 


15582 






















24.84 


2.34 


16984 






















24.85 


4.24 


15854 






















24.86 


0.89 


7907 








18 


10.77 


0.92 


3263 








24.80 


3.75 


2486 










10.78 


0.88 


6944 








24.81 


3.69 


315 










10.79 


0.88 


5944 








24.82 


3.62 


15585 










10.80 


1.8 


4564 








24.83 


3.68 


18409 






















24.84 


4.23 


18376 






















24.85 


3.62 


3861 






















24.86 


2.70 


4175 








20 


10.77 


2.8 


2850 














25.42 


1.46 


5533 




10.78 


2.9 


6445 














25.44 


0.90 


11482 




10.79 


2.8 


8192 














25.46 


0.64 


15773 




10.80 


4.4 


5806 














25.48 


2.31 


15951 






















25.50 


1.44 


15117 


22 


10.77 


2.1 


3621 


11.06 


3.1 


5472 
















10.78 


2.1 


6737 


11.08 


2.9 


7019 
















10.79 


2.1 


10700 






















10.80 


3.4 


5256 




















24 


10.77 


0.12 


1452 


11.06 


1.08 


5959 
















10.78 


2.9 


9063 


11.08 


1.1 


9821 
















10.79 


2.8 


13614 






















10.80 


2.8 


4475 




















28 








11.06 


1.5 


6753 






















11.08 


1.4 


11935 














30 








11.06 


1.2 


2308 






















11.08 


1.2 


15251 















TABLE VL The statistics (in millions of composite sweeps) used in the re-weighting analysis for (3c in SU(A'c) gauge theory. 
Also quoted is the integrated auto-correlation time. Studies at Nt = 10 and 12 have been carried out with statistics of 4 x 10® 
composite sweeps, where Tint varied between 500 and 1000 sweeps. 
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_i 

Ns = 16 

10.40 0.502033(46 
10.46 0.520680(25 
10.48 0.525249(12 
10.50 0.529260(25 
10.52 0.532786(15 
10.54 0.535885(9) 
10.60 0.543796(10 
10.70 0.554088(7 
10.76 0.559352(8 
10.78 0.561000(4 
10.79 

10.80 0.562606(7 
10.82 0.564167(7 
10.84 0.565700(8 
10.90 0.570122(5 
11.00 0.576972(5 
11.02 0.578291(5 
11.04 0.579577(4 
11.06 0.580859(5 
11.08 0.582117(5 
11.10 0.583361(5 
11.20 0.589361(4 
11.30 0.595051(3 
11.50 0.605665(5 
11.60 0.610655(6 



SU(4) 

= 18 
0.501991(29) 
0.520728(23) 
0.525254(16) 
0.529268(17) 
0.532749(20) 
0.535903(16) 
0.543778(12) 
0.554097(10) 

0.561005(10) 
0.561805(8) 

0.562605(5) 
0.564166(8) 
0.565689(5) 
0.570123(4) 
0.576979(5) 



0.583369(3) 
0.589357(3) 
0.595043(4) 
0.605652(4) 
0.610643(4) 



iV. = 24 
0.502136(18 
0.520673(16 
0.525259(7) 
0.529274(15 
0.532773(11 
0.535903(12 
0.543790(7 
0.554092(4 
0.559341(5 
0.561000(5 
0.561799(4 
0.562602(3 
0.564169(4 
0.565702(3 
0.570120(3 
0.576984(3 
0.578286(4 
0.579580(3 
0.580851(4 
0.582109(4 
0.583363(3 
0.589360(2 
0.595051(2 
0.605650(2 
0.610630(3 



SU(6) 



N, = 16 



jP) 
N, =20 



iV, =24 



24 


60 0.537686(30) 






24 


70 0.542241(24) 






24 


80 


0.546406(9) 






24 


82 0.547189(10) 






24 


84 


0.547981(6) 






24 


86 


0.548764(9) 






24 


88 


0.549532(8) 






24 


90 


0.550292(7) 






24 


92 


0.551042(6) 






25 


00 


0.553961(3) 






25 


10 


0.557459(5) 






25 


20 


0.560805(4) 






25 


30 


0.564026(4) 


0.564012(4) 


0.564015(3 


25 


40 


0.567126(4) 


0.567119(2) 


0.567120(3 


25 


42 




0.567727(3) 


0.567727(2 


25 


44 




0.568332(3) 


0.568329(3 


25 


46 




0.568933(3) 


0.568931(3 


25 


48 




0.569522(3) 


0.569527(3 


25 


50 


0.570130(5) 


0.570123(4) 


0.570119(3 


25 


52 




0.570709(3) 


0.570708(2 


25 


54 




0.571293(4) 


0.571294(2 


25 


56 




0.571880(3) 


0.571877(3 


25 


60 


0.573043(3) 


0.573031(3) 


0.573031(2 


25 


70 


0.575876(4) 






25 


80 


0.578631(3) 


0.578620(3) 


0.578616(3 


25 


90 


0.581321(4) 






26 


00 


0.583937(5) 


0.583925(4) 


0.583921(4 


26 


10 


0.586506(6) 


0.586487(3) 


0.586485(3 



TABLE VII; Plaquette expectation values for SU(4) and SU(6). The numbers in brackets denote errors on the least significant 
digits. Note that the volume dependence is negliglible. 



SU(8) 


SU(IO) 


/3 {P) 


(3 (P) 


44.50 0.542133(9) 
44.80 0.548772(6) 
45.00 0.552860(7) 
45.20 0.556747(5) 

45.50 0.562252(4) 
46.00 0.570739(4) 


68.00 0.376818(6) 
70.00 0.542542(8) 
71.00 0.554874(8) 
72.00 0.567207(8) 
74.00 0.586693(5) 
76.00 0.603381(7) 



TABLE VIII: Plaquette expectation values for SU(8) and SU(IO), measured on 16* lattices. The numbers in brackets are errors 
on the least significant digit. 



16 



SU(4) 


SU(6) 




/3 




Nt = 


6 


Nt = 


8 




P 


Nt=6 


Nt=8 


10 


7000 


0.9106(9) 





6963 


1) 


24 


6000 


0.8841(4) 


0.6714(8) 


10 


7200 


0.9311(9) 





7119 


1) 


24 


7000 


0.9335(5) 


0.7089(9) 


10 


7400 





9513 


10) 





7274 


1) 


24 


8000 


0.9818(5) 


0.7456(9) 


10 


7600 





9716 


10) 





7429 


1) 


24 


8200 


0.9913(5) 


0.7529(9) 


10 


7700 





9817 


10) 





7507 


1) 


24 


8400 


1.0010(5) 


0.7602(9) 


10 


7800 





9920 


10) 





7585 


1) 


24 


8500 


1.0059(5) 


0.7640(9) 


10 


7900 


1 


0020 


10) 





7662 


1) 


24 


8600 


1.0107(5) 


0.7675(9) 


10 


8000 


1 


0122 


10) 





7740 


1) 


24 


8800 


1.0204(5) 


0.7749(9) 


10 


8200 


1 


0326 


10) 





7895 


1) 


24 


9000 


1.0301(5) 


0.7823(9) 


10 


8400 


1 


0530 


11) 





8052 


1) 


24 


9200 


1.0397(5) 


0.7896(10) 


10 


9000 


1 


1151 


11) 





8526 


1) 


25 


0000 


1.0786(5) 


0.8191(10) 


11 


0000 


1 


2215 


12) 





9340 


1) 


25 


1000 


1.1277(6) 


0.8564(10) 


11 


0200 


1 


2433 


13) 





9506 


1) 


25 


2000 


1.1775(6) 


0.8943(11) 


11 


0400 


1 


2654 


13) 





9675 


1) 


25 


3000 


1.2282(6) 


0.9327(11) 


11 


0600 


1 


2876 


13) 





9845 


1) 


25 


3400 


1.2487(6) 


0.9483(11) 


11 


0800 


1 


3101 


13) 


1 


0017 


1) 


25 


3800 


1.2695(6) 


0.9641(12) 


11 


1000 


1 


3331 


13) 


1 


0193 


1) 


25 


4000 


1.2799(6) 


0.9720(12) 


11 


1200 


1 


3559 


14) 


1 


0368 


1) 


25 


4200 


1.2904(6) 


0.9800(12) 


11 


1400 


1 


3791 


14) 


1 


0545 


1) 


25 


4400 


1.3008(6) 


0.9879(12) 


11 


1600 


1 


4028 


14) 


1 


0726 


1) 


25 


4600 


1.3114(7) 


0.9960(12) 


11 


1800 


1 


4265 


14) 


1 


0908 


1) 


25 


4800 


1.3221(7) 


1.0040(12) 


11 


2000 


1 


4507 


15) 


1 


1092 


1) 


25 


5000 


1.3327(7) 


1.0121(12) 


11 


2200 


1 


4748 


15) 


1 


1276 


1) 


25 


5200 


1.3434(7) 


1.0202(12) 


11 


2400 


1 


4996 


15) 


1 


1466 


1) 


25 


5400 


1.3542(7) 


1.0284(12) 


11 


2600 


1 


5243 


15) 


1 


1655 


1) 


25 


5600 


1.3650(7) 


1.0366(13) 


11 


2800 


1 


5496 


16) 


1 


1849 


1) 


25 


6000 


1.3868(7) 


1.0532(13) 


11 


3000 


1 


5755 


16) 


1 


2046 


1) 


25 


8000 


1.4989(7) 


1.1383(14) 


11 


3200 


1 


6010 


16) 


1 


2242 


1) 


25 


9000 


1.5573(8) 


1.1827(14) 


11 


3400 


1 


6271 


16) 


1 


2441 


1) 


26 


0000 


1.6169(8) 


1.2280(15) 


11 


3600 


1 


6537 


17) 


1 


2645 


1) 


26 


2000 


1.7415(9) 


1.3226(16) 


11 


3800 


1 


6804 


17) 


1 


2849 


1) 


26 


4000 


1.8727(9) 


1.4222(17) 


11 


4000 


1 


7077 


17) 


1 


3058 


1) 


26 


6000 2.0119(10) 


1.5279(18) 


11 


4200 


1 


7351 


17) 


1 


3267 


1) 


27 


0000 2.3149(11) 


1.7581(21) 


11 


4400 


1 


7629 


18) 


1 


3479 


1) 


27 


5000 2.7474(14) 


2.0865(25) 


11 


4600 


1 


7911 


18) 


1 


3696 


1) 










11 


4800 


1 


8197 


18) 


1 


3914 


1) 










11 


5000 


1 


8486 


19) 


1 


4135 


1) 










11 


6000 


1 


9988 


20) 


1 


5283 


1) 










11 


7000 


2 


1590 


22) 


1 


6508 


1) 










11 


8000 


2 


3300 


23) 


1 


7816 


1) 










11 


9000 


2 


5126 


25) 


1 


9212 


1) 










12 


0000 


2 


7076 


27) 


2 


0703 


1) 










12 


2000 


3 


1390 


32) 


2 


4001 


1) 










12 


4000 


3 


6320 


37) 


2 


7771 


1) 











TABLE IX: T/Tc scales for SU(4) and SU(6) gauge theories in the V-scheme, for Nt = 6 and 
the uncertainty in the determination of /3c. 



8. The errors are dominated by 



